function GAD_Mid_cdf_plot_rm(Data_Summary)
%% WT saline data
G18_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F4_18') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G18_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G18_saline_RM = RM;
        G18_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F4_18') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G18_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G18_saline_RM = cat(1,G18_saline_RM,RM);
        G18_saline_Shuffle = cat(1,G18_saline_Shuffle,RM_Shuffle);
    end
end
[G18sf,G18sx,G18slow,G18shigh] = ecdf(G18_saline_RM,'Function','cdf','Alpha',0.05);
[G18ssf,G18ssx,G18sslow,G18sshigh] = ecdf(G18_saline_Shuffle,'Function','cdf','Alpha',0.05);

G28_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G28_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_saline_RM = RM;
        G28_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G28_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_saline_RM = cat(1,G28_saline_RM,RM);
        G28_saline_Shuffle = cat(1,G28_saline_Shuffle,RM_Shuffle);
    end
end
[G28sf,G28sx,G28slow,G28shigh] = ecdf(G28_saline_RM,'Function','cdf','Alpha',0.05);
[G28ssf,G28ssx,G28sslow,G28sshigh] = ecdf(G28_saline_Shuffle,'Function','cdf','Alpha',0.05);

G31_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G31_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_saline_RM = RM;
        G31_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G31_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_saline_RM = cat(1,G31_saline_RM,RM);
        G31_saline_Shuffle = cat(1,G31_saline_Shuffle,RM_Shuffle);
    end
end
[G31sf,G31sx,G31slow,G31shigh] = ecdf(G31_saline_RM,'Function','cdf','Alpha',0.05);
[G31ssf,G31ssx,G31sslow,G31sshigh] = ecdf(G31_saline_Shuffle,'Function','cdf','Alpha',0.05);

G34_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G34_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_saline_RM = RM;
        G34_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G34_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_saline_RM = cat(1,G34_saline_RM,RM);
        G34_saline_Shuffle = cat(1,G34_saline_Shuffle,RM_Shuffle);
    end
end
[G34sf,G34sx,G34slow,G34shigh] = ecdf(G34_saline_RM,'Function','cdf','Alpha',0.05);
[G34ssf,G34ssx,G34sslow,G34sshigh] = ecdf(G34_saline_Shuffle,'Function','cdf','Alpha',0.05);

%% Organize data for WT saline

% get WT saline CDF mean
A = cat(1,G18sf,G28sf,G31sf,G34sf);
WT_saline_F = sort(A);

% get WT saline CDF x-steps
B = cat(1,G18sx,G28sx,G31sx,G34sx);
WT_saline_X = sort(B);

% get WT saline CDF high bound mean
C = cat(1,G18shigh,G28shigh,G31shigh,G34shigh);
WT_saline_H = sort(C);

% get WT saline CDF low bound mean
D = cat(1,G18slow,G28slow,G31slow,G34slow);
WT_saline_L = sort(D);

%% MUT saline data
G17_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F4_17') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G17_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G17_saline_RM = RM;
        G17_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F4_17') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G17_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G17_saline_RM = cat(1,G17_saline_RM,RM);
        G17_saline_Shuffle = cat(1,G17_saline_Shuffle,RM_Shuffle);
    end
end
[G17sf,G17sx,G17slow,G17shigh] = ecdf(G17_saline_RM,'Function','cdf','Alpha',0.05);
[G17ssf,G17ssx,G17sslow,G17sshigh] = ecdf(G17_saline_Shuffle,'Function','cdf','Alpha',0.05);

G21_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G21_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_saline_RM = RM;
        G21_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G21_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_saline_RM = cat(1,G21_saline_RM,RM);
        G21_saline_Shuffle = cat(1,G21_saline_Shuffle,RM_Shuffle);
    end
end
[G21sf,G21sx,G21slow,G21shigh] = ecdf(G21_saline_RM,'Function','cdf','Alpha',0.05);
[G21ssf,G21ssx,G21sslow,G21sshigh] = ecdf(G21_saline_Shuffle,'Function','cdf','Alpha',0.05);

G45_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G45_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_saline_RM = RM;
        G45_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G45_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_saline_RM = cat(1,G45_saline_RM,RM);
        G45_saline_Shuffle = cat(1,G45_saline_Shuffle,RM_Shuffle);
    end
end
[G45sf,G45sx,G45slow,G45shigh] = ecdf(G45_saline_RM,'Function','cdf','Alpha',0.05);
[G45ssf,G45ssx,G45sslow,G45sshigh] = ecdf(G45_saline_Shuffle,'Function','cdf','Alpha',0.05);

G53_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G53_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_saline_RM = RM;
        G53_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G53_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_saline_RM = cat(1,G53_saline_RM,RM);
        G53_saline_Shuffle = cat(1,G53_saline_Shuffle,RM_Shuffle);
    end
end
[G53sf,G53sx,G53slow,G53shigh] = ecdf(G53_saline_RM,'Function','cdf','Alpha',0.05);
[G53ssf,G53ssx,G53sslow,G53sshigh] = ecdf(G53_saline_Shuffle,'Function','cdf','Alpha',0.05);


%% Organize data for MUT saline

% get MUT saline CDF mean
A = cat(1,G17sf,G21sf,G45sf,G53sf);
MUT_saline_F = sort(A);

% get MUT saline CDF x-steps
B = cat(1,G17sx,G21sx,G45sx,G53sx);
MUT_saline_X = sort(B);

% get MUT saline CDF high bound mean
C = cat(1,G17shigh,G21shigh,G45shigh,G53shigh);
MUT_saline_H = sort(C);

% get MUT saline CDF low bound mean
D = cat(1,G17slow,G21slow,G45slow,G53slow);
MUT_saline_L = sort(D);

%% WT etom 7mg/kg data
G18_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F4_18') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G18_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G18_etom7_RM = RM;
        G18_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F4_18') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G18_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G18_etom7_RM = cat(1,G18_etom7_RM,RM);
        G18_etom7_Shuffle = cat(1,G18_etom7_Shuffle,RM_Shuffle);
    end
end
[G18e7f,G18e7x,G18e7low,G18e7high] = ecdf(G18_etom7_RM,'Function','cdf','Alpha',0.05);
[G18se7f,G18se7x,G18se7low,G18se7high] = ecdf(G18_etom7_Shuffle,'Function','cdf','Alpha',0.05);

G28_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G28_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_etom7_RM = RM;
        G28_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G28_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_etom7_RM = cat(1,G28_etom7_RM,RM);
        G28_etom7_Shuffle = cat(1,G28_etom7_Shuffle,RM_Shuffle);
    end
end
[G28e7f,G28e7x,G28e7low,G28e7high] = ecdf(G28_etom7_RM,'Function','cdf','Alpha',0.05);
[G28se7f,G28se7x,G28se7low,G28se7high] = ecdf(G28_etom7_Shuffle,'Function','cdf','Alpha',0.05);

G31_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G31_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_etom7_RM = RM;
        G31_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G31_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_etom7_RM = cat(1,G31_etom7_RM,RM);
        G31_etom7_Shuffle = cat(1,G31_etom7_Shuffle,RM_Shuffle);
    end
end
[G31e7f,G31e7x,G31e7low,G31e7high] = ecdf(G31_etom7_RM,'Function','cdf','Alpha',0.05);
[G31se7f,G31se7x,G31se7low,G31se7high] = ecdf(G31_etom7_Shuffle,'Function','cdf','Alpha',0.05);

G34_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G34_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_etom7_RM = RM;
        G34_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G34_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_etom7_RM = cat(1,G34_etom7_RM,RM);
        G34_etom7_Shuffle = cat(1,G34_etom7_Shuffle,RM_Shuffle);
    end
end
[G34e7f,G34e7x,G34e7low,G34e7high] = ecdf(G34_etom7_RM,'Function','cdf','Alpha',0.05);
[G34se7f,G34se7x,G34se7low,G34se7high] = ecdf(G34_etom7_Shuffle,'Function','cdf','Alpha',0.05);

%% Organize data for WT Etom 7mg/kg

% get WT etom7 CDF mean
A = cat(1,G18e7f,G28e7f,G31e7f,G34e7f);
WT_etom7_F = sort(A);

% get WT etom7 CDF x-steps
B = cat(1,G18e7x,G28e7x,G31e7x,G34e7x);
WT_etom7_X = sort(B);

% get WT etom7 CDF high bound mean
C = cat(1,G18e7high,G28e7high,G31e7high,G34e7high);
WT_etom7_H = sort(C);

% get WT etom7 CDF low bound mean
D = cat(1,G18e7low,G28e7low,G31e7low,G34e7low);
WT_etom7_L = sort(D);

%% MUT etom 7mg/kg data
G17_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F4_17') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G17_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G17_etom7_RM = RM;
        G17_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F4_17') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G17_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G17_etom7_RM = cat(1,G17_etom7_RM,RM);
        G17_etom7_Shuffle = cat(1,G17_etom7_Shuffle,RM_Shuffle);
    end
end
[G17e7f,G17e7x,G17e7low,G17e7high] = ecdf(G17_etom7_RM,'Function','cdf','Alpha',0.05);
[G17se7f,G17se7x,G17se7low,G17se7high] = ecdf(G17_etom7_Shuffle,'Function','cdf','Alpha',0.05);

G21_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G21_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_etom7_RM = RM;
        G21_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G21_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_etom7_RM = cat(1,G21_etom7_RM,RM);
        G21_etom7_Shuffle = cat(1,G21_etom7_Shuffle,RM_Shuffle);
    end
end
[G21e7f,G21e7x,G21e7low,G21e7high] = ecdf(G21_etom7_RM,'Function','cdf','Alpha',0.05);
[G21se7f,G21se7x,G21se7low,G21se7high] = ecdf(G21_etom7_Shuffle,'Function','cdf','Alpha',0.05);

G45_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G45_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_etom7_RM = RM;
        G45_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G45_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_etom7_RM = cat(1,G45_etom7_RM,RM);
        G45_etom7_Shuffle = cat(1,G45_etom7_Shuffle,RM_Shuffle);
    end
end
[G45e7f,G45e7x,G45e7low,G45e7high] = ecdf(G45_etom7_RM,'Function','cdf','Alpha',0.05);
[G45se7f,G45se7x,G45se7low,G45se7high] = ecdf(G45_etom7_Shuffle,'Function','cdf','Alpha',0.05);

G53_etom7_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')  && isempty(G53_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_etom7_RM = RM;
        G53_etom7_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'etomidate 7mg/kg')   && ~isempty(G53_etom7_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_etom7_RM = cat(1,G53_etom7_RM,RM);
        G53_etom7_Shuffle = cat(1,G53_etom7_Shuffle,RM_Shuffle);
    end
end
[G53e7f,G53e7x,G53e7low,G53e7high] = ecdf(G53_etom7_RM,'Function','cdf','Alpha',0.05);
[G53se7f,G53se7x,G53se7low,G53se7high] = ecdf(G53_etom7_Shuffle,'Function','cdf','Alpha',0.05);


%% Organize data for MUT Etom 7mg/kg

% get MUT etom7 CDF mean
A = cat(1,G17e7f,G21e7f,G45e7f,G53e7f);
MUT_etom7_F = sort(A);

% get MUT etom7 CDF x-steps
B = cat(1,G17e7x,G21e7x,G45e7x,G53e7x);
MUT_etom7_X = sort(B);

% get MUT etom7 CDF high bound mean
C = cat(1,G17e7high,G21e7high,G45e7high,G53e7high);
MUT_etom7_H = sort(C);

% get MUT etom7 CDF low bound mean
D = cat(1,G17e7low,G21e7low,G45e7low,G53e7low);
MUT_etom7_L = sort(D);

%% Finally, put shuffle data together

% get shuffle CDF mean
A = cat(1,G17ssf,G21ssf,G45ssf,G53ssf,G18ssf,G28ssf,G31ssf,G34ssf,G17se7f,G18se7f,G21se7f,G28se7f,G31se7f,G34se7f,G45se7f,G53se7f);
Shuffle_F = sort(A);
Shuffle_F = cat(1,0,Shuffle_F,1);

% get shuffle CDF x-steps
B = cat(1,G17ssx,G21ssx,G45ssx,G53ssx,G18ssx,G28ssx,G31ssx,G34ssx,G17se7x,G18se7x,G21se7x,G28se7x,G31se7x,G34se7x,G45se7x,G53se7x);
Shuffle_X = sort(B);
Shuffle_X = cat(1,-1,Shuffle_X,1);

% get shuffle CDF high bound mean
C = cat(1,G17sshigh,G21sshigh,G45sshigh,G53sshigh,G18sshigh,G28sshigh,G31sshigh,G34sshigh,G17se7high,G18se7high,G21se7high,G28se7high,G31se7high,G34se7high,G45se7high,G53se7high);
Shuffle_H = sort(C);

% get shuffle CDF low bound mean
D = cat(1,G17sslow,G21sslow,G45sslow,G53sslow,G18sslow,G28sslow,G31sslow,G34sslow,G17se7low,G18se7low,G21se7low,G28se7low,G31se7low,G34se7low,G45se7low,G53se7low);
Shuffle_L = sort(D);

%% plot all data

figure
subplot(1,4,1)
plot(WT_saline_X,WT_saline_F,'Color','#1750AC','LineWidth',3);
hold on
% % plot(WT_saline_X,WT_saline_L,'--k','LineWidth',1);
% % hold on
% % plot(WT_saline_X,WT_saline_H,'--k','LineWidth',1);
% % hold on
plot(MUT_saline_X,MUT_saline_F,'Color','#F53BD6','LineWidth',3);
hold on
% % plot(MUT_saline_X,MUT_saline_L,'--k','LineWidth',1);
% % hold on
% % plot(MUT_saline_X,MUT_saline_H,'--k','LineWidth',1);
% % hold on
plot(Shuffle_X,Shuffle_F,'-k','LineWidth',3);
hold on
% % plot(Shuffle_X,Shuffle_L,'--k','LineWidth',1);
% % hold on
% % plot(Shuffle_X,Shuffle_H,'--k','LineWidth',1);
% % hold on
% % legend('p-WT Saline','','','\alpha5-i-KO Saline','','','Shuffle','95% Confidence Interval','','Location','southeast');
xlabel('RM Correlation Value')
ylabel('Cumulative Fraction of Cells')
xlim([-0.6,1])
ylim([0,1])
set(gca,'FontSize',10)
set(gca, 'visible', 'off')
hold off

subplot(1,4,2)
plot(WT_saline_X,WT_saline_F,'Color','#1750AC','LineWidth',3);
hold on
% % plot(WT_saline_X,WT_saline_L,'--k','LineWidth',1);
% % hold on
% % plot(WT_saline_X,WT_saline_H,'--k','LineWidth',1);
% % hold on
plot(WT_etom7_X,WT_etom7_F,'Color','#73B9EE' ,'LineWidth',3);
hold on
% % plot(WT_etom7_X,WT_etom7_L,'--k','LineWidth',1);
% % hold on
% % plot(WT_etom7_X,WT_etom7_H,'--k','LineWidth',1);
% % hold on
plot(Shuffle_X,Shuffle_F,'-k','LineWidth',3);
hold on
% % plot(Shuffle_X,Shuffle_L,'--k','LineWidth',1);
% % hold on
% % plot(Shuffle_X,Shuffle_H,'--k','LineWidth',1);
% % hold on
% % legend('p-WT Saline','','','p-WT Etom 7mg/kg','','','Location','southeast');
xlabel('RM Correlation Value')
ylabel('Cumulative Fraction of Cells')
xlim([-0.6,1])
ylim([0,1])
set(gca,'FontSize',10)
set(gca, 'visible', 'off')
hold off

subplot(1,4,3)
plot(MUT_saline_X,MUT_saline_F,'Color','#F53BD6','LineWidth',3);
hold on
% % plot(MUT_saline_X,MUT_saline_L,'--k','LineWidth',1);
% % hold on
% % plot(MUT_saline_X,MUT_saline_H,'--k','LineWidth',1);
% % hold on
plot(MUT_etom7_X,MUT_etom7_F,'Color','#FA9CEA' ,'LineWidth',3);
hold on
% % plot(MUT_etom7_X,MUT_etom7_L,'--k','LineWidth',1);
% % hold on
% % plot(MUT_etom7_X,MUT_etom7_H,'--k','LineWidth',1);
% % hold on
plot(Shuffle_X,Shuffle_F,'-k','LineWidth',3);
hold on
% % plot(Shuffle_X,Shuffle_L,'--k','LineWidth',1);
% % hold on
% % plot(Shuffle_X,Shuffle_H,'--k','LineWidth',1);
% % hold on
% % legend('\alpha5-i-KO Saline','','','\alpha5-i-KO Etom 7mg/kg','','','Location','southeast');
xlabel('RM Correlation Value')
ylabel('Cumulative Fraction of Cells')
xlim([-0.6,1])
ylim([0,1])
set(gca,'FontSize',10)
set(gca, 'visible', 'off')
hold off

subplot(1,4,4)
plot(WT_etom7_X,WT_etom7_F,'Color','#73B9EE' ,'LineWidth',3);
hold on
% % plot(WT_etom7_X,WT_etom7_L,'--k','LineWidth',1);
% % hold on
% % plot(WT_etom7_X,WT_etom7_H,'--k','LineWidth',1);
% % hold on
plot(MUT_etom7_X,MUT_etom7_F,'Color','#FA9CEA' ,'LineWidth',3);
hold on
% % plot(MUT_etom7_X,MUT_etom7_L,'--k','LineWidth',1);
% % hold on
% % plot(MUT_etom7_X,MUT_etom7_H,'--k','LineWidth',1);
% % hold on
plot(Shuffle_X,Shuffle_F,'-k','LineWidth',3);
hold on
% % plot(Shuffle_X,Shuffle_L,'--k','LineWidth',1);
% % hold on
% % plot(Shuffle_X,Shuffle_H,'--k','LineWidth',1);
% % hold on
% % legend('p-WT Etom 7mg/kg','','','\alpha5-i-KO Etom 7mg/kg','','','Location','southeast');
xlabel('RM Correlation Value')
ylabel('Cumulative Fraction of Cells')
xlim([-0.6,1])
ylim([0,1])
set(gca,'FontSize',10)
set(gca, 'visible', 'off')
hold off


end